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VLF  Workshop:  24-25  January  1285 


Marine  Physical  Laboratory 
Scripps  Institution  of  Oceanography 
San  Diego,  CA 


I.  Introduction 

The  focus  of  this  classified  workshop  was  on  the  frequency  region  1-30  Hi.  The  program  included 
presentations  by  the  participants  followed  by  a  discussion  covering  areas  of  potential  Navy  interest  and 
future  research  objectives.  A  summary  of  the  presentations  and  discussion  is  contained  in  this  report. 


II.  Participants 

The  following  is  a  roster  of  those  people  who  participated  in  the  VLF  Workshop.  Full  addresses 
and  telephone  numbers  are  provided  in  Appendix  A. 
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Mr.  James  McEachern 

Naval  Postgraduate  School 

Dr.  Rudolph  Nichols 

APL/JH 

Dr.  Julius  Bowen 
Mr.  Joe  Lombardo 
Mr.  Ross  Rottier 

ARL/UT 


MIT/WHOI 

Dr.  Arthur  Baggeroer 
Mr.  John  Ewing 

SIO 

Dr.  LeRoy  Dorman 
Dr.  John  Orcutt 
Dr.  Peter  Worcester 

MPL/SIO 

Dr.  Wm.  Hodgkiss 
Dr.  Victor  Anderson 
Dr.  Fred  Fisher 
Mr.  Lee  Culver 
Mr.  Chris  Nickles 
Mr.  Greg  Edmonds 

PSI 

Dr.  Marshall  Bradley 

BBN 

Mr.  Mike  Sullivan 
Mr.  Dale  Green 
Mr.  Alan  Ma 

Polar  Research,  Inc. 

Dr.  James  Wilson 


Dr.  Robert  Koch 


ID.  Presentation* 

The  presentations  covered  the  following  broad  areas:  (l)  ambient  noise,  (2)  propagation,  and  (3) 
system  implications.  Summaries  of  the  presentations  (i.e.  abstract  and  copies  of  the  vu- graphs)  are 
contained  either  in  Appendix  B  of  this  report  or  in  a  companion  classified  MPL  technical  report.  The 
presentations  were  given  in  the  following  order. 

Ambient  Noise 

1.  W.S.  Hodgkiss  and  V.C.  Anderson,  N  Infrasonic  Acoustic  Ambient  Noise" 

2.  Art  Baggeroer  (for  Ira  Dyer),  "Arctic  Ambient  Noise" 

3.  Marshall  Bradley,  "  VLF  Ambient  Noise" 

4.  Jim  Wilson,  "Active  Ice  Ridge  Ambient  Noise  Levels" 

5.  John  Orcutt,  "  Noise  Characteristics  and  Wave  Propagation  at  Low  Frequencies" 

Propagation 

1.  John  Orcutt,  "Noise  Characteristics  and  Wave  Propagation  at  Low  Frequencies"  (continued) 

2.  LeRoy  Dorman,  "  The  Seafloor  Sound  Channel  and  Stoneley  Waves  as  Observed  Using 
Seafloor  Explosions  and  Ocean  Bottom  Seismographs" 

3.  Jim  Wilson,  "  Buck/ Wilson  Deep  Water  Arctic  Transmission  Loss  Model" 

4.  Art  Baggeroer,  "  Energy  Partitioning  for  Long  Range,  Low  Frequency  Propagation  in  the 
Arctic  Ocean" 

5.  John  Ewing,  "VLF  Propagation  in  the  Oceanic  Crust  and  Upper  Mantle" 

6.  Art  Baggeroer,  "  On  the  Relative  Amplitudes  between  Primary  and  Multiple  Signals  from 
Seismic  Refractions  in  Oceanic  Crust" 

7.  Ralph  Baer,  "A  Three-Dimensional  Model  for  Bathymetric  Scattering  with  VLF 
Applications" 

8.  Tom  Tunnell,  "Modeling  Support  for  Cape  Fear  VLF  Exercise" 

System  Implications 

1.  Jim  McEachern,  "  NADC  Infrasonic  Sonobuoy  Program  -  A  Brief  Summary"  * 

2.  Victor  C.  Anderson,  "  Comments  on  Nearfield  VLF  Surveillance" 

3.  Julius  Bowen,  "Effect  of  Submarine  Speed  on  Comparative  Detectability  of  Blade  and 
Machinery  Tones"  * 

4.  Roes  Rottier,  "  Transmission  Loss  to  a  Seismometer  in  Oceanic  Basement" 

5.  Richard  Koch,  "VLF  Propagation  Modeling  to  an  Ocean  Borehole  Receiver" 

8.  Joe  Lombardo,  "Applications  of  Ocean  Seismometers  to  ASW"  ** 

7.  Wm.  Carey,  "VLF  Borehole  Surveillance  Program  Requirements"  * 

8.  Mike  Morrison,  "Arctic  Remote  Sensing"  * 

9.  Randy  Moore,  "Soviet  Submarine  Source  Levels"  * 

10.  Mike  Sullivan,  "  Implications  of  VLF  Clutter  on  Performance  Prediction"  * 

*  All  or  a  subset  of  the  vu-graphs  are  contained  in  a  companion  classified  MPL  technical  report. 


**  Vu-graphs  not  available. 


IV.  Summary 

The  last  half  day  of  the  workshop  was  devoted  to  reviewing  and  summarizing  what  had  been  said 
during  the  course  of  the  presentations.  In  general,  the  discussion  was  a  blend  of  the  three  broad  areas 
covered  in  the  workshop:  (1)  ambient  noise,  (2)  propagation,  and  (3)  system  implications. 

A.  Ambient  Noise 

Based  on  available  data  and  current  models,  the  characteristics  of  ambient  noise  in  the  VLF 
region  appear  to  be  dominated  by  a  microseism  peak  in  the  region  of  0.1  Hz  with  shipping  noise 
being  the  dominant  contributor  above  3-4  Hz.  More  needs  to  be  understood  regarding  ambient 
noise  levels  generated  due  to  atmospheric  pressure  fluctuations,  surface  wave  motion,  and 
turbulence.  Future  ambient  noise  measurement  programs  should  acquire  a  substantial  amount  of 
environmental  data  to  aid  in  isolating  the  various  contributing  mechanisms  (meterological,  water 
column,  and  bottom  properties).  Additional  array  studies  are  needed  to  characterize  the 
frequency/wavenumber  content  of  VLF  noise  in  order  to  better  understand  where  the  noise  is 
coming  from  and  what  propagation  paths  are  important. 

B.  Propagation 

At  VLF  frequencies,  bottom  interaction  is  a  significant  aspect  of  the  propagation  problem. 
There  is  a  strong  need  to  have  available  range-dependent  propagation  models  for  use  in 
investigating  the  significant  pathways  of  energy  transfer.  For  example,  future  models  based  on 
modal  formalisms  must  have  the  capability  of  properly  accounting  for  the  coupling  of  energy 
between  modes  in  order  to  explain  existing  observations.  High  quality  array  studies  at  these  low 
frequencies  are  needed  to  accurately  measure  amplitudes,  phase  velocities,  and  angles  of  approach 
of  propagating  energy.  Many  complex  wave  propagation  phenomena  can  be  adequately  modeled 
with  algorithms  which  simply  take  into  account  all  multiples.  Nevertheless,  strong  evidence  exists 
to  support  an  important  role  for  scattering  at  these  low  frequencies  although  the  incorporation  of 
a  useful  theory  in  the  advanced  propagation  formalisms  remains  an  unsolved  problem. 

C.  System  Implications 

A  better  understanding  of  the  origin  and  magnitude  of  VLF  ambient  noise  coupled  with 
measured/predicted  source  levels  and  new  propagation  models  will  aid  in  ascertaining  the 
usefulness  of  the  VLF  band  for  surveillance  purposes.  The  trade-oils  between  the  various  types  of 
sensors  (geophone/hydrophone)  and  their  placement  (water  column,  bottom  mounted,  or  buried) 
need  to  be  better  understood.  Sensor  design  to  reduce  flow-induced  self-noise  and  the  decoupling 
of  tethered  sensors  from  the  contaminating  effects  of  mechanical  linkages  (strumming)  are  areas  of 
continuing  interest.  Over  what  aperture  the  coherent  processing  of  data  from  an  array  of  VLF 
sensors  is  viable  (coherence  length)  and  whether  the  array  should  be  vertical  or  horizontal  are 
important  system-related  issues. 
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Appendix  B 

This  appendix  contains  summaries  (i.e.  abstracts  and  vu-graphs)  of  most  of  the  presentations.  As 
indicated  earlier,  summaries  of  a  few  of  the  presentations  are  contained  in  a  companion  classified  MPL 
technical  report. 
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Fig.  7 

This  spectrogram  shows,  among  other  things,  a  77  Hz  ice  event 
lasting  -for  more  than  20  minutes.  It  is  influenced  by  the 
increased  ice  stress  caused  by  the  helicopter  landing  at  time  7 
minutes,  and  by  the  decreased  ice  stress  caused  by  the  underwater 
air  gun  at  16  minutes. 


FIG.  1 


Composite  of  ambient  noise  typical  of  the  central  Arctic  under 
moderate  noise  conditions.  Measurements  made  at  a  depth  of  90  m 
below  the  ice.  The  frequency  character i sti c  below  1  Hz  is  not 
yet  explained,  but  hypotheses  include  nonlinear  surface  wave 
noise  from  the  open  ocean,  pseudosound  from  the  planetary 
boundary  layer  turbulence  interacting  with  the  hydrophone,  and 
radiation  from  the  seismic  noise  in  the  earth's  crust  and  mantle. 
Peaks  from  1-10  Hz  are  caused  by  hydrophone  cable  strum,  which  is 
quite  variable  even  under  similar  cable  and  current  conditions. 
The  broad  peak  centered  on  about  15  Hz  is  most  likely  caused  by  a 
large  number  of  independent  ice  cracking  events. 

Fig.  2 

This  shows  in  clearer  form  the  broad  peak  centered  at  about  15  Hz 
mentioned  in  Fig.  1,  A  60  Hz  tone  due  to  line  frequency  pick  up 
has  been  eliminated. 

Fig.  3 

A  time  series  of  noise  (pressure  magnitude  in  the  10-20  Hz  band) 
is  shown  for  a  10  day  period  in  the  central  Arctic.  (Noise 
magnitude  is  sampled  once  per  hour  in  a  10  minute  window.) 

Fig.  4 

Inferred  ice  stress  far  the  same  10  day  period  as  in  Fig.  3. 

This  ice  stress  is  the  horizontal  normal  stress  in  the  ice 
determined  as  a  consequence  of  force  balance  among  wind  loads, 
current  loads,  coriolis  forces,  and  pressure  gradient  forces. 

This  inferred  internal  ice  stress  is  followed  closely  by  the  time 
series  of  Fig.  3.  (Ice  stress  is  sampled  once  per  hour  in  a  10 
minute  window.) 

Fig.  5 

Cross-correlation  of  the  time  series  of  Figs.  3  and  4,  indicating 
that  horizontal  internal  ice  stress  is  a  good  predictor  of 
acoustic  noise.  Other  data  show  that  the  noise  is  created  by  the 
large  number  of  independent  ice  cracking  events,  with  the  most 
likely  mechanism  being  flexural  failure  of  the  ice.  These  ice 
cracking  events  are  individually  on  the  order  of  a  few  100  msecs 
in  duration. 

Fig.  6 

A  pathological  observation  in  the  marginal  ice  zone  and  often 
seen  as  well  in  the  central  Arctic.  This  pathology  represents  a 
very  long  time  narrow  band  event  lasting  at  times  for  minutes  to 
hours.  They  are  hypothesized  to  be  caused  by  slip-stick 
oscillations  of  adjacent  floes  in  shear  deformation. 
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Fig.  3.  Surface  echo  data:  July  27,  1982.  (Edge  detection  algorithm 

input) 


Fig.  4.  Error  time  series:  July  27,  1982.  (Edge  detection  algorithm 
output) 


this  observation.  The  resulting  windowed-edge  detection 
time  series  of  sea  surface  range  corresponding  roughly  to 
8-33  h  in  Fig.  3  is  shown  in  Fig.  S  (1.S  min/sample).  A  1024- 
point  FFT  of  that  time  series  is  provided  in  Fig.  6.  Period  (in¬ 
verse  frequency)  is  displayed  along  the  horizontal  axis  for 
convenience.  The  broad  band  of  power  from  approximately 
120  to  30  min  in  period  corresponds  to  internal  wave  activity. 
Note  that  the  spectrum  has  dropped  off  to  being  essentially 
flat  beyond  IS  min  in  period.  Thus  localizing  the  array  ele¬ 
ments  at  least  as  often  as  every  7.5  min  would  capture  the 
significant  dynamics  of  array  motion. 

IV.  CONCLUSIONS 

Surface  echo  data  from  a  recent  sea  test  of  the  prototype 
element  of  a  freely  drifting  infrasonic  measurement  system  has 
been  used  to  illustrate  the  type  of  pulse  processing  which  will 
be  implemented  as  a  first  step  in  the  localization  procedure. 
For  the  surface  echo  data  set  processed,  an  arrival  time  error 
of  1 .3-ms  rms  was  indicated  by  the  edge  detection  algorithm. 
Since  this  represents  an  uncertainty  of  approximately  0.1  X  at 
10  Hz,  we  can  expect  our  overall  rms  positioning  eiror  for 
the  full  array  to  be  less  than  0.1  X  [12] .  Such  random  posi¬ 
tioning  errors  lead  to  a  maximum  array  signal  gain  degradation 
oflessthan2dB  (13). 


Fig.  5.  Surface  range  time  series:  July  27, 1982.  (Edge  detection  algo¬ 
rithm  output) 

. . . . . . . . . . 1 
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Fig.  6.  FFT  of  surface  range  time  series:  July  27, 1982.  (Edge  detec¬ 
tion  algorithm  output) 
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depicted  in  Fig.  2.  This  8-kHz  system  has  a  maximum  range 
of  approximately  IS  fan  and  a  nominal  resolution  of  ±0.75  m 
based  on  a  1-ms  sampling  interval  at  the  output  of  the  pulse 
detection  circuit.  A  different  buoy  pings  every  1 5  min  until 
each  Swallow  float  has  transmitted.  Hie  sequence  then  re¬ 
peats.  Since  each  buoy  maintains  an  internal  time  base  (clock), 
deciphering  which  buoy  transmitted  any  given  received  pulse 
is  not  difficult.  Synchronization  of  the  clocks  is  possible  due 
to  the  occurrence  of  reciprocal  transmission  paths  between 
buoys.  All  floats  (including  the  transmitter  listening  to  its 
own  sea  surface  reflection)  internally  record  arrival  time  of  the 
transmission  on  magnetic  tape.  These  tapes  then  will  be  read* 
at  the  completion  of  a  deployment  and  the  arrival  times  at 
all  buoys  from  each  ping  will  be  combined  algorithmically  to 
yield  the  relative  three-dimensional  positions  of  the  array  ele¬ 
ments  as  a  function  of  time.  An  absolute  coordinate  system 
can  be  established  with  the  addition  of  external  receivers 
and/or  transmitters  of  known  position  (eg.,  a  ship  positioned 
via  satellite  navigation  and  a  bottom-tethered  buoy  at  a  known 
location,  as  shown  in  Fig.  2). 

III.  TIME-OF-ARRIVAL  ESTIMATION 

Within  each  buoy,  the  arrival  times  of  the  received  pings 
must  be  estimated  as  a  first  step  in  the  localization  process. 
The  pulse-detection  circuit  is  of  a  fairly  common  design.  It 
consists  of  two  filters  in  cascade.  Both  pass  the  localization 
pulse  which  is  centered  at  8  kHz  and  is  6  ms  in  duration.  The 
first  filter  is  rather  broad  band  (2-kHz  bandwidth  centered  at 
10  JS  kHz).  The  output  of  this  filter  is  hard  limited  yielding  a 
constant  power  time  series.  The  second  filter  is  narrow  band 
(0.2-kHz  bandwidth  centered  at  8  kHz).  The  power  spectrum 
of  ambient  noise  passed  on  to  the  second  filter  predominately 
lies  outside  its  bandwidth  and  thus  yields  a  small  output.  In 
contrast,  the  power  spectrum  of  an  8-kHz  pulse  received  on 
the  skirt  of  the  first  filter  mostly  lies  inside  the  bandwidth  of 
the  second  filter  and  thus  yields  a  large  output.  The  envelope- 
detected  output  of  the  second  filter  is  compared  against  a 
threshold  and  generates  a  pulse  present-absent  decision.  The 
decision  output  [1,  0]  of  the  sonar  is  sampled  at  a  I-kHz 
rata.  That  bit  stream  is  not  recorded  directly,  but  is  first  oper¬ 
ated  on  by  data -compression  logic.  The  bit  stream  is  inspected 
in  blocks  of  8  bits.  If  a  1  (pulse  present)  is  observed,  the  block 


of  8  bits  is  recorded  as  one  byte  along  with  its  associated  time 
(2  bytes).  Otherwise,  nothing  is  recorded.  Up  to  85  3-byte 
groups  can  be  recorded  for  each  1 .5-min  period. 

Both  ocean-ambient  noise  and  gross-data  gaps  make  the  raw 
sonar  data  unusable  as  a  continuous  record  of  pulse  time  of 
arrival.  Accurate  arrival  time  measurements  are  required  by 
array-element  localization  and  tracking  algorithms.  Thus  it  is 
necessary  to  preprocess  the  raw  detector  output  to  provide 
usable  inputs  to  these  algorithms.  We  have  had  several  oppor¬ 
tunities  to  test  our  prototype  self-contained  Swallow  float  at 
sea.  These  tests  have  provided  experimental  data  sets  repre¬ 
sentative  of  those  which  will  be  recorded  by  the  acoustic 
localization  system.  Using  one  of  these  data  sets,  we  have 
investigated  candidate  algorithms  for  the  preprocessing  opera¬ 
tion. 

As  mentioned  previously,  each  buoy  not  only  listens  for 
transmissions  from  the  other  elements  of  the  array,  but  for  the 
sea  surface  reflection  of  its  own  pulse  as  well.  Raw  sonar 
surface  echo  data  from  a  recent  sea  trip  is  shown  in  Fig. 
3.  Time  delay  after  transmission  is  displayed  positive  down¬ 
ward  in  terms  of  range.  The  closest  edge  of  the  dark  band 
represents  the  shortest  path  to  the  surface  from  the  Swallow 
float.  Hence,  the  problem  of  distinguishing  the  surface  echo 
from  other  returns  and  noise  becomes  a  problem  of  detecting 
this  edge. 

Several  edge  detection  algorithms  have  been  tried.  The 
following  approach,  which  makes  use  of  an  adaptive  linear 
predictor  directing  the  search  for  an  edge  in  a  narrow  range 
window,  has  been  the  most  successful. 

First,  range  is  predicted  using  a  16-weight  one-step  linear 
predictor  incorporating  the  least  mean-squares  (LMS)  adaptive 
algorithm  [11].  Then,  the  detector  proceeds  through  four 
modes  until  a  range  is  decided. 

1)  If  no  data  record  is  present,  the  predicted  range  is  used. 

2)  If  a  data  record  is  present,  a  search  is  made  for  an  edge 
in  a  small  window  plus  and  minus  two  standard  devia¬ 
tions  of  the  error  (i.e.,  rms  error)  from  the  predicted 
range. 

3)  If  no  edge  is  found  in  the  small  window,  a  larger  window 
representing  the  maximum  physically  possible  move¬ 
ment  is  searched. 

4)  If  the  large  window  search  fails,  the  predicted  range  is 
used. 

After  a  range  is  determined,  the  edge  detector  calculates  the 
error  and  the  adaptive  predictor  is  iterated. 

The  use  of  a  small  window  is  based  on  the  assumption  that 
the  range  will  not  deviate  much  from  the  predicted  value.  By 
using  a  multiple  of  the  rms  error  for  the  width  of  the  small 
window,  the  window  is  allowed  to  adapt  to  increasing  or 
decreasing  noise  in  the  data.  The  large  window  allows  any 
physically  reasonable  movement  to  be  accepted,  even  though 
it  might  not  be  predictable. 

When  used  with  a  small  window  of  two  standard  deviations, 
a  large  window  of  10  ms,  and  the  July  27,  1982  data  set,  the 
edge  detector  remained  in  mode  2  most  of  the  time,  indicating 
the  predictor  was  predicting  the  range  well.  As  shown  in 
Fig.  4,  the  small  error  signal  (rms  error  of  1 .3  ms)  confirmed 
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Abstract  -Initial  testing  of  the  prototype  element  of  a  freely  drifting 
infrasonic  sensor  array  is  described.  Tbe  intent  of  this  measurement 
system  is  to  gather  wide  aperture  data  sets  which  will  be  need  both  to 
characterize  ambient  noise  in  the  region  1-10  Hz  and  to  assess  the 
gains  possible  from  beam  forming  utilizing  a  collection  of  eery  low 
frequency  (VLF)  sensors.  Coherent  processing  (beam  forming)  of  the 
Infrasonic  sensor  data  is  made  possible  by  relative  position  measure¬ 
ments  derived  from  mutual  acooatic  interrogation  of  the  elements  at  a 
higher  frequency.  Surface  echo  data  from  a  recent  sea  test  of  the 
prototype  buoy  are  used  to  Oloatrate  the  type  of  poise  processing 
which  will  be  implemented  as  a  first  step  in  the  localization  procedure. 

I.  INTRODUCTION 

NE  PROJECT  at  the  Marine  Physical  Laboratory  (MPL)  is 
the  development  of  an  array  to  study  infrasonic  acoustic 
background  noise  in  the  ocean.  In  general,  self-noise  limits  an 
accurate  assessment  of  the  true  nature  of  this  ambient  acoustic 
noise  field.  Cable  strumming,  flow  noise,  and  the  local  pressure 
and  velocity  fields  due  to  turbulence  are  all  difficulties  associ¬ 
ated  with  measurements  in  the  region  1-10  Hz  when  using  an* 
array  which  requires  some  form  of  tethering  or  mechanical 
linkage  [l]-[9].  Even  for  deep  sensors  near  or  on  the  sea 
floor,  the  flow  noise  associated  with  very  slow  currents  can 
significantly  contaminate  the  local  pressure  or  particle  velocity 
field. 

The  approach  taken  in  this  project  is  to  develop  an  auton¬ 
omous  buoy  capable  of  recording  the  components  of  particle 
velocity  in  the  1-10- Hz  band  and  which  is  also  able  to  both 
generate  and  receive  high-frequency  acoustic-positioning 
signals.  The  buoy  is  a  neutrally  buoyant  Swallow  float  which  is 
stable  at  midwater  depths  and,  drifting  freely,  is  not  subject  to 
any  flow  disturbance.  The  acoustic-positioning  capability  will 
be  used  to  monitor  the  relative  geometry  of  a  set  of  these 
buoys  in  a  freely  drifting  infrasonic  array,  so  that  data  from 
the  individual  elements  can  be  combined  coherently  off-line 
by  use  of  the  MPL  dynamic  beamformer  [  10] . 

A  prototype  buoy  has  been  designed,  built,  and  deployed 
at  sea.  The  performance  of  the  acoustic-positioning  system  in 
this  prototype  buoy,  in  measuring  the  time  delay  to  the  first 
surface  echo,  provides  insight  into  the  localization  accuracy 
that  can  be  expected  for  a  multi-element  array  system.  Satis¬ 
factory  performance  of  the  acoustic-positioning  system  is  a 
key  factor  in  the  success  of  the  project. 

Manuscript  nceivsd  December  29,  1982;  revised  April  20,  1983. 
This  work  wu  supported  by  the  Office  of  Naval  Research  under  Con¬ 
tract  N00014-80-C-0220.  Contribution  of  the  Scripps  Institution  of 
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II.  THE  INFRASONIC  SENSOR  BUOY  AND  ARRAY 

The  element  of  the  Swallow  float  array  is  a  spherical  glass 
Swallow  float  which  contains  three  geophones  that  measure 
the  three  components  of  particle  velocity  in  the  1-10-Hz 
band:  a  compass  for  buoy  heading;  an  acoustic  transponder  for 
localization;  a  solid-state  memory-data  buffer;  a  digital -tape 
data  recorder;  and  an  acoustically  actuated  ballast  release.  The 
general  configuration  of  hardware  in  our  prototype  unit  is 
shown  in  Fig.  1 .  The  floats  are  neutrally  buoyant  and  can  be 
ballasted  for  a  desired  depth.  Limited  by  tape  recorder  capac¬ 
ity,  the  maximum  submergence  period  of  the  floats  is  on  the 
order  of  60  h.  Each  buoy  continuously  fills  an  8-kB  buffer 
memory,  then  periodically  writes  this  buffer  out  to  tape.  The 
information  contained  in  the  buffer  includes  compass  heading, 
acoustic  localization  pulse  arrival  tL.ies,  and  three  channels  of 
geophone  data  sampled  at  25  Hz.  Currently,  the  buffer  is 
written  out  to  tape  in  a  burst  mode  every  1 .5  min.  During  this 
1-s  procedure,  no  data  are  sampled. 

Several  of  these  elements  can  be  deployed  in  a  random 
configuration  to  create  a  spatially  distributed  array.  Their 
relative  three-dimensional  position  as  a  function  of  time  are 
obtained  via  an  acoustic  pulse  mutual  interrogation  system  as 
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Potential  Sources  of  Measure  lent  Contamination 

o  Cable  Strumming 
o  Flow  Noise 

o  Turbulence-Related  Pressure  Fluctuations 


Motivating  Questions 

What  is  the  temporal  and  spatial  structure  of  ambient  ocean 
noise  in  the  1-10  Hz  frequency  region? 

How  do  the  results  from  different  sensor  types  compare  (e.g. 
VLF  sonobuoys,  bottom  tethered  hydrophones,  and  bottom 
mounted  hydrophones)? 
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Lack  of  quality  VLF  data 

No  generally  accepted  theory  for  the  non 
shipping  noise  component  near  10  Hz 

No  routine  predictive  capability 

Great  potential  for  simple  hydrophones 
being  self  noise  limited 


Hydrophones  mounted  on  cables  are  likely 
to  be  contaminated  by  strum  noise 
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HOTEL  II  3  5  500  377  x  30  x  25  2  30  000 


SOURCES  OF  VLF  AMBIENT  NOISE 


Seismic  Sources 

Microseisms 

Shipping 

Surface  Waves 

Wave  Turbulence  Interaction  with  Surface  Waves 
Wind  Turbulence 
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Turbulent  Pressure  Fluctuations 


Current 


Seismic  interference  in  the  VLF  band 

•  Seismic  profiling  -  Channel  1 

•  T-Phase  -  Channel  9,  10,  11 


Seismic  interference  in  the  VLF  band 

•  Seismic  profiling  -  Channel  1,  5 

•  Signals  -  Channel  4,  6,  7 


sets  momeirer  on 
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SPECTRUM  LEVEL  (dB/»  Pt/Hi'*) 


VLF  Sound  Radiation  by  Merchant  Ships 
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Received  levels  at  100  miles  due  to  the  freighter 
ADOLF  LEONHARDT  (22,000  tons,  10,600  bhp)  moving 
at  15  knots  (from  Wittenborn,  1976)  . 
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Figure*  2  -  P^n  View  of  the  Experiment 
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CONCLUSIONS 


In  deep  water  areas  shipping  is  usually 
dominant  in  the  10-20  Hz  band 

Below  about  4  Hz  the  sharp  increase  in  noise 
levels  appears  to  be  caused  by  microseisms 

Difficult  to  assign  generation  mechanisms  for 
the  hydrodynamic  component  of  VLF  noise 

No  convincing  comparisons  between  measurement 
and  theory 

Noise  levels  in  the  4-10  Hz  band  are  comparable 
to  levels  in  the  10-20  Hz  band 


No  predictive  capability 


Nonacoustic  Noise  Interference  in  Measurements 
of  Infrasonic  Ambient  Noise  (Strassberg,  1979) 

Mechanisms 

Hydrophone  response  to  temperature 
inhomogeneities 

Turbulence  generated  by  flow  past 
hydrophone 

•  Interaction  between  hydrophone  and 
turbulence 


Wind  Component  of  VLF  Noise 


Dependence  of  the  non- shipping  component  of  VLF  ambient 
noise  upon  wind  speed  (from  Wittenborn,  1976) . 
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Figure  8  -  Receiver  Depth  Dependence  Predicted  by  FACT 
for  Various  Source  Depths. 


Figure  10  -  Comparison  of  Estimated  Noise  Spectra  and 
Measured  Noise  Spectra  In  the  Beaufort  Sea 


John  Orcutt 
Rick  Adair 
Tom  Sereno 
Peter  Shearer 

•  Seafloor  and  subseafloor  noise 

0  Signal  comparisons  on  seafloor  and  subseaf loor  Inertial 
and  hydrophone  Instruments  at  various  phase  velocities 

0  Propagation  of  setsmlc/acoustic  signals  over  large 
ranges  <Pn.  Sn,  T) 


8  Subseafloor  noise  levels  as  much  as  20  db  lower  than 
seafloor  noise  levels. 

0  Seafloor  noise  levels  consistent  with  many  previous 
acoustics  measurements. 

0  Seafloor  noise  at  low  frequencies,  and  high  frequencies  In 
the  presence  of  shipping.  Is  coherent  over  baselines  of  at  least 
700  m. 

8  Seafloor  noise  dominated,  to  4-5  Hz,  by  Interaction  of 
surface  gravity  waves. 

0  Seafloor  noise  above  5  Hz  likely  controlled  by  shipping. 

0  Horizontal  noise  at  seafloor  gener  ally  greater  than 
verttcal  noise.  The  levels  are  equal  at  subseafioor  receivers. 

0  Signal  levels  at  subseafloor  sites  are  comparable  In  size 
to  seafloor  sites.  Signal  to  noise  ratios  Increase  with  burial  of 
sensor. 

•  Pn,  Sn,  T  phases  can  be  synthesized  by  wave  propagation 

In  a  laterally  homogeneous  medium  with  no  inter-lithosphere 
waveguides. 

0  Slow  decay  of  Pn,  Sn,  T  phases  with  range  depends  upon 

the  high  Q  ocean  and  not  upon  a  high  Q  within  the  crust  and 
lithosphere. 

0  Scattering  plays  a  minor  role  In  the  formation  of  long  Pn, 
Sn,  T  wavetralns. 

0  The  character  of  the  Pn,  Sn,  T  wavetraln  Is  controlled  by 
the  thickness  of  the  water  and  sediment  columns. 
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FIGURE  1.  Velocity  and  density  model.  The  left  side  shows 
entire  model  while  on  the  right  the  depth  scale  is  expanded 
to  show  detail  in  the  crust.  Asterisks  indicate  source  depths 
The  model  is  from  Gettrust  and  Frazer  (1981). 
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FIGURE  3.  Complete  synthectic  seismograms  for  a  thrust  fault 
source  at  14.0  km  depth  at  an  epicentral  range  of  1000  km 
a)  Vertical  component,  b)  horizonfal  component  rotated  25 
degrees  clockwise  from  the  radial  direction.  Velocities  of 
8.1,  4.7,  2.6  and  1.6  km/sec  are  indicated. 
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FIGURE  4.  Synthetic  Pn  wavetrain.  Contains  the  first 
100  seconds  of  the  records  in  figure  3.  a)  Vertical 
component  and  b)  horizontal  component. 
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'IGURE  5.  Record  section  of  the  vertical  normal  stress,  crzz0 
at  14.0  km  depth.  Lines  are  drawn  with  an  inverse  slope 
of  8.17  km/sec  separated  by  the  two-way  travel  time  in  the 
water  column,  5.33  sec. 
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FIGURE  6.  Power  spectrum  of  a  reverberating  wavelet,  a) 
WavetrairLjeomprised  of  a  reverberating  wavelet  with  a 
spacing  between  multiples  of  T=0.46  seconds  and  atten¬ 
uation  factor,  a  =  1.0.  b)  Power  spectrum  of  the  wavetrain 
in  (a)  superimposed  on  the  spectrum  of  a  single  wavelet 
(dashed).  The  peaks  are  equally  spaced  with  df=1/T=2.17 
hz  and  are  offset  by  an  amount  df/2. 
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URE  9.  Spectral  characteristics  of  sediment  and  water 
iverberations.  a)  Wavelet,  w(t);  b)  wavelet  spectrum,  |W(f)|2 
)  Sediment  reverberation  withl_zxs=3.0  and  Ts=0.086  sec; 

/  sediment  reverberation  spectrum,  e)  Water  reverberation 
ith  aw=0.05  and  Tw=8.0  sec;  f)  first  4.0  hz  of  the  water 
jverberation  spectrum,  g)  Full  wavetrain,  x(t);  h)  wavetrain 
pectrum,  |X(f)|2. 
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FIGURE  10.  Full  wavetrain  spectrum,  |X(f)|2.  a)  First  3.0  hz 
of  the  spectrum  with  vertical  lines  at  fn=(2n  +  1  )dfw/2, 
(n=0,1 ,2,...).  b)  Complete  spectrum  with  vertical  lines  at 
fn=(2n+1  )dfs/2,  (n=0,1,2,..  ). 


11.  Power  spectra  of  the  synthetic  Pn  wavetrains  in  figure 
3  vertical  component  spectrum  is  shown  in  (a)-(c)  and  the 
mtal  component  spectrum  in  (d)-(f).  Figures  (c)  and  (f) 
le  complete  spectra.  Figures  (b)  and  (e)  were  obtained  by 
ation  of  a  running  mean  filter  to  the  complete  spectra, 
d  lines  indicate  the  predicted  frequencies  associated  with  .; 
ent  reverberation,  Iona  bars  for  S-waves  and  short  bars  to 
/es.  Figures  (a)  and  (d)  are  the  first  3.0  hz  of  the  spectra 
Jashed  lines  at  fn=(2n  +  1  )df/2  (n=0, 1 ,2,...)  where  df  is 
id  by  equation  (11). 
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FIGURE  15.  OBS  data,  vertical  component,  a)  OBS  Juan,  ever 
14,  b)  OBS  Karen,  event  24,  c)  OBS  Suzy,  event  25,  d)  91 
Suzy,  event  100  and  e)  OBS  Suzy,  event  130. 


FIGURE  16.  OBS  data,  hydrophone,  a)  OBS  Suzy,  event  35,  b) 
OBS  Karen,  event  24,  c)  OBS  Suzy,  event  25,  d)  OBS  Suzy, 
event  100  and  e)  OBS  Suzy,  event  130. 


FIGURE  17.  The  effect  of  the  instrument  on  the  synthetic  Pn 
wavetrain.  (a)  Synthetic  Pn  wavetrain  of  figure  4a  with  no 
instrument  and  (b)  with  the  vertical  component  OBS  responst. 
mcluded.  (c)  OBS  Suzy,  event  100  vertical  component,  (d) 
synthetic  Pn  wavetrain  with  hydrophone  response  included  and 
(e)  OBS  Suzy,  event  100;  hydrophone. 
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FIGURE  14.  Schematic  representation  of  the  isolation  of  a 
low  frequency  spectral  peak  through  the  constructive 
interference  of  different  layer  reverberations.  The  figures 
represent  idealized  spectra,  a)  Wavelet  spectrum,  b) 
P-wave  sediment  reverberation  spectrum,  c)  S-wave  sed 
iment  reverberation  spectrum  and  d)  wavetrain  spectrum. 
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CHANNEL  VERT 
MAX:  1.68E+01 
Y-TIC:0.SE+01 
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NO  FILTER 
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(PASS  BAND) 
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CHANNEL  VERT 
MAX:  9-27E+00 
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1.0-  4  0  Hz 
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CHANNEL  VERT 
MAX:  1.00E+01 
Y-T1C:0.2E+01 
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FIGURE  19.  Bandpassed  versions  of  the  seismoc-;>ir 
15(c).  The  seismogram  in  (a)  is  unfilterea  and  t 
bands  in  (b)  — (d)  are  0  25  to  1.0  hz,  1.0  to  4.C 
4.0  to  16.0  hz  respectively. 
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FIGURE  20.  Power  spectrum  of  the  seismogram  in  figure 
15(b).  Dashed  lines  indicate  predicted  positions  of  the 
spectral  peaks  associated  with  P-wave  reverberations  i 
the  sediments. 
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FIGURE  21.  First  3.0  hz  of  the  power  spectrum  of  the 
seismogram  in  figure  16(e).  Vertical  lines  indicate 
the  predicted  positions  of  the  spectral  peaks  corres¬ 
ponding  to  water  reverberations. 


OBSERVATIONS  OF  STONELEY  WAVES  USING  SEA  FLOOR  SHOTS  AND  RECEIVERS 
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Anthony  E.  Schreiner 

The  spatial  and  temporal  correlation  properties  of  background  noise  are 
primarily  controlled  by  two  things:  1)  the  excitation  and  2)  the  waveguides 
in  which  the  noise  travels.  We  focus  our  attention  here  on  the  waveguide  at 
the  sea  floor.  We  have  excited  and  recorded  Stoneley  waves  along  the  sea 
floor  in  order  to  extract  the  physical  properties  of  the  top  few  meters  of  the 
sea  floor.  Our  location  was  the  site  of  DSDP  hole  #469,  in  3800  m  of  water  at 
the  base  of  the  Patton  escarpment  west  of  San  Diego.  We  used  2.5  kg  charges 
of  TNT  fired  by  electronic  timers.  Five  of  the  eight  shots  attempted  per¬ 
formed  satisfactorily.  Three  records,  at  ranges  between  500  and  1500  m  showed 
Stoneley  waves,  while  shots  at  greater  ranges  did  not.  Fitting  the  observed 
group  velocities  (or  slownesses)  to  a  simple  model  by  trial  and  error  yielded 
a  model  with  a  linear  shear  velocity  gradient  of  3.5s  and  shear  velocity  at 
the  sea  floor  of  about  50  ms 

We  also  measured  the  response  of  the  sea  floor  to  excitation  by  the  in¬ 
strument  and  find  that  data  consistent  with  a  damped  system  with  a  resonant 
frequency  of  18  Hz  and  a  Q  of  5. 

Both  the  physical  property  work  and  the  resonance  measurement  bear  on 
the  problem  of  equalization  of  OBS  responses  so  that  they  may  be  used  pro¬ 
ductively  in  arrays. 


Multichannel  seismic-reflection  profile  AA'  showing  proposed  position  of  Site  469  (S.  P.  Lee).  (The  location  of  this  profile  is  shown  on  Fig.  I.) 


THUMPER  CRU] 


-120.58  -120.58  -120.54 


-120.58  -120.56  -120.54 


LONG' 


he  least  squares  solution 
triangles  are  OBSs  wh i 1 


TL  TL0(f>R)  +  lOlogR  +  ATL( 


zs>zr'f)  +  a(f,o) 


Short  Range  Spreading 
"57. 50+. 10(f-10)+aQ(f )R  10Hz<f<20Hz 

58. 50+ . 50(f-20)' 5+aQ(f )R  20Hz<f<200Hz 

.65.20+. 01(f-2 00)+  a  (f)R  200Hz<f<S00Hz - 
TLQ(f)  =  min  Long  Range  Spreading 

"62 . 00+  .30(f-10)  10Hz<f<20Hz 

65 . 00+ . 65 ( f-20 ) ‘ 6  20Hz<f<200Hz 

,79.65+.01(f-200)  200Hz<f<300hz  . 


ATL  =  ATL  =-10  log  sin2(27Tsin0zf /c) 
r  s 


ATL  = 


"  . 5zf /c- . 72  10Hz^f<50Hz 

( .9-.008f )zf/c-(1.296-.01152f )  50Hzsf<100  Hz 

(  .lzf/c-.144)(^|^2  100Hz<f^800  Hz 


a(  f  ,  a)  -  CQ(f)0n(f) 


n(f) 


2f2+550 
f  2+50 


-9 

f4 . 32 

10Hz<i<20Hz 

2.35x10 

2 . 40xl0-6 

f2.00 

20Hz<f<100Hz 

2. 40xl0'3 

f0.500 

100Hz<f<200Hz 

-2 

3.40x10 

200Hz<f<  BOOHz 

J  -cr 


TL  REPORT 


PAGE  2 


quantitative  difference  in  the  values  of  these  functions  are  found, 

because  the  assumptions  about  the  nature  of  ice  scattering  is  very 

different  in  different  models.  The  Buck/Wilson  ice  scattering  function 

sVMyrf-  j 

derived  in  this  paper  is  divided  into  two  parts  for  -shot-  and  long 
ranges  and  given  a  physical  interpretati on. 

The  development  of  the  TL  model  and  analysis  of  the  data  uncovered 
a  possible  "megaphone"  or  slope  enhanced  propagation  effect  on  the 
eastern  slopes  of  the  Northwind  Ridge.  It  is  emphasized  that  the 
megaphone  effect  is  very  site  dependent  and  no  implications  are  made 
for  the  e>si stance  of  the  megaphone  effect  over  other  sloping  areas. 
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ABSTRACT 


A  semi -empirical  Arctic  transmission  Loss  <TL)  model  is  developed 
based  on  TL  data  collected  between  ice  camps  and/or  aircraft  in  the 
deep  Arctic  basins  since  1970.  The  model  is  limited  to  -frequencies 
between  10Hz  and  QOOH2 ,  to  source  and  receiver  depths  above  BOO  ft., 
and  to  minimum  water  depths  between  source  and  receiver  of  100m  or 
greater.  The  model  estimates  average  TL  and  no  attempt  is  made  to 
account  for  spread  in  TL  at  fixed  ranges. 

The  Buck/Wilson  Arctic  TL  Model  consists  of  four  components:  the 
source/receiver  depth  dependence  of  TL;  cyclindrical  spreading; 
spherical  spreading  plus  low  (<10°)  angle  and  bottom  interaction  loss; 
and  long  range,  ice  scattering  attenuation.  Since  measured  TL  data 
connot  be  uniquely  divided  into  these  four  areas,  care  is  taken  to  use 
basic  TL  properties  known  to  be  true  in  the  open  ocean.  The  depth 
dependence  of  TL  is  modeled  using  standard  surface/image  interference 
functions  and  the  data  show  that  the  pressure  release  surface  is  the 
sea/ice  interface  and  not  the  air/ice  interface.  The  parabolic 
equation  <PE)  and  the^SEPS  Transmission  Loss  (ASTRAL)  model  estimates 
of  spreading  loss  are  used  as  a  guide  in  deriving  spreading  losr 
functions.  The  "residual"  TL,  defined  as  the  measured  TL  minus  the 
model  TL  depth  and  spreading  functions,  is  used  to  derive  the 
attenuation  due  to  ice  scattering  loss  as  a  function  of  frequency  and 

f' 

ice  roughness  (standard  deviation).  The  comparison  of  the  drived  ice 
scattering  attenuation  with  other  model  ice  scattering  attenuation  or 
reflection  functions  will  be  adressed  in  a  'future  paper.  A  largo 


SCUFFS  INSTITUTION  OF  OCEiNOGRAHY  OCEAN  BOTTOM  SEISMOGRAPHS 


RECORDER: 

TAPE: 

TAPE  FORMAT: 

CAPACITY: 

MEMORY: 

DYNAMIC  RANGE: 

NODES: 

SCHEDULE  PROGRAMMING: 
STATUS: 

INTERNAL  TIMING: 
POSITIONING: 

POWER: 

DIMENSIONS: 

WEIGHT: 

REFERENCE: 


10M  r*el-to  reel  1/4"  tap* 

7200  f**t. 

4  tracks  digital  aerial  biphase.  System  status 
and  time  distributed  with  the  data.  3  seismometer 
channels.  1  hydrophone. 

24  hoars  recording. 

5.2  minutes  for  4  components  at 
128  samples  per  second  per  channel. 

102  DB  (12  bit  ADC  with  64  x  high  speed 
gain  ranging). 

Programmed  and/or  triggered 

FI ezible 

Message  containing  noise  level,  tape  nsage,  battery 
voltages  transmitted  to  surface  on  acoustic  cosmiand. 

TCXO  to  60  days 

Acoustic  transponder 

About  400  MW 

About  1  m*  footprint.  22  "i.d.  sphere  with  6" 
center  ring. 

About  550  pounds  in  air,  with  ballast  tripod. 

Moore,  R.D. .  L.M.  Dorman,  C-Y  Huan,  and  D.L.  Berliner. 
An  Ocean  bottom,  microprocessor  based  seismometer. 
Marine  Geophysical  Researches,  4,  451-477,  1981. 
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The  relationship  between  group  velocity  of  a  mod*  and  ampli¬ 
tude  resulting  -from  a  point  fore*  excitation.  It  shows  that 
th#  largest  spectral  amplitude  occurs  at  the  minimum  of  the 
group  velocity  of  the  mode. 
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Figure  63  Relative  Energy  of  Deep  and  Shallow  RSR  Rays  for  4  Source 
Depths  (Average  for  all  4  Bottom  Categories) 
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TITLE  FREQUENCY  SOURCE  RECEIVER 
ASTRAL  10.0  HZ  300-0  FT  800-0  FT 
ASTRAL  50-0  HZ  300-0  FT  800-0  FT 
ASTRAL  200-0  HZ  300-0  FT  800-0  FT 
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VLF  PROPAGATION  #1 


ENERGY  PARTITIONING  FOR  LONG  RANGE,  LOW  FREQUENCY  PROPAGATION  IN 

THE  ARCTIC  OCEAN 


The  partitioning  o-f  energy  ■from  explosive  signals 
propagating  in  the  Arctic  Ocean  'for  ranges  up  to  350  km  was 
determined  as  a  -function  o-f  -frequency  using  a  1  km  x  I  km  array 
during  the  FRAM  ice  station  experiments.  The  energy  was  resolved 
into  three  classes:  i)  dispersive,  low  order  modes  (1-3)  trapped 
in  the  upward  re-fracting  Arctic  sound  channel;  ii)  multiple  RSR 
paths;  and  iii)  and  extensive  set  o-f  bottom  re-fracted  signals  on 
a  reversed  branch  o-f  the  travel  time  curve.  Estimating  the 
energy  of  each  component  as  a  -function  o-f  travel  time  and  ray 
parameter  indicates  that  below  20  Hz,  the  bottom  interacting 
components  contain  well  over  hal-f  the  total  observed  energy  when 
the  bathymetry  between  the  source  and  receiver  is  smooth. 
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Sonogram  analyslo  of  channel  1.  Contour  Interval  ■  3  dB. 

Solid  curve  n  initial  model'  )  Dotted  curve  -  corrected  model. 
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VLF  Propagation  in  the  Oceanic  Crust  and  Upper  Mantle 


Abstract 


Studies  of  the  characteristics  of  VLF  (1-5  Hz)  propagation  at  ranges  greater  than  30  km  show 
considerable  coherency  both  in  primary,  intracrustal,  and  water  column  multiple  refraction  energy.  The 
data  utilized  were  generated  by  large  explosive  charges  and  recorded  by  Ocean  Bottom  Hydrophones  on 
the  East  Pacific  Rise  during  the  ROSE  project,  along  with  some  comparable  data  obtained  from  the 
Mid  Atlantic  Ridge.  The  water-crust-upper-mantle  wave  guide  behaves  in  an  orderly  and  predictable 
fashion  that  accounts  for  the  distribution  of  energy  among  the  primary  and  multiple  paths.  Preliminary 
tests  using  corsscorrelation  between  nearby  seismograms  suggests  that  for  a  number  of  range  windows, 
in  which  multipath  influence  does  not  seriously  degrade  the  data,  coherency  is  sufficiently  high  to  permit 
signal  summation  of  array  elements  distributed  over  an  area  with  dimensions  as  large  as  10  km. 


A'%  > 


Locations  of  receivers  and  large  shots  for  ROSE  data  study.  OBH  seismogram 
sections  recorded  at  1,  6,  7,  8,  and  8',  indicated  by  arrows  are  shown  in 
figures.  Dashed  line  shows  location  of  the  East  Pacific  Rise  axis. 


PU  £L.  Velocity-depth  model  determined  from  shots  and  Instruments  along 
Line  4,  from  Bratt  and  Purdy.  In  press.  Dashed  line  below  6.7  km 
Indicates  velocity  structure  appropriate  for  upper  mantle  for 
propagation  perpendicular  to  ridge  axis. 


Figure^:  Composite  seismogram  section  of  shots  on  Line  1  toward 
the  west  from  Instruments  6,  7,  and  8.  Solid  lines  are 
drawn  through  primary  Pn  and  multiple  path  Pn  arrival 
times;  dashed  lines  connect  the  primary  PmP  and  multiple 
path  PmP  arrivals.  Examples  of  propagation  paths  are 
shown  for  Pn,  WPn,  2WPn,  and  W2Pn  arrivals  in  Figure  3. 
Note  the  shift  of  peak  energy  to. higher-order  multiples 
at  ranges  of  50  and  75  km. 


Pn+2WPn 


Propagation  paths  for  several  primary  and  multiple  phases  observed 
in  the  data.  Phase  designations  are  the  same  as  those  in  Figure  2 
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^‘5  $.  Seismogram  section  recorded  at  OBH  8'  on  Line  4.  Shots  are  north 
of  the  Instrument.  Mote  the  absence  of  Pn  energy  and  the  dominance 
of  primary  PmP  and  multi  ole  PmP  phases. 
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BATHYMETRIC  CONTOUR  CHAR 
MULTI-BEAM  DATA^J 


U.S.  NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  DC  20375 


DEPARTMENT  OF  THE  NAVY 
NAVAL  RESEARCH  LABORATORY 
WASHINGTON,  O.C.  20375 
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A  Three-Dimensional  Model  for  Bathymetric  Scattering 
with  VLF  Applications 


R.N.  Baer,  J.S.  Perkins,  E.B.  Wright,  and  B.B.  Adams 


We  have  developed  a  computer  model  to  simulate  the  performance  of 
acoustic  arrays  receiving  energy  scattered  from  rough  ocean  bathymetry. 

The  model  has  been  applied  using  a  30  Hz  source  In  the  vicinity  of  a 
seamount  where  bathymetric  statistics  are  available.  Predictions  have 
been  obtained  for  a  long  horizontal  array  and  a  short  vertical  array.  We 
discuss  the  degrading  effects  of  the  rough  bathymetry  (Including  the  high 
sidelobe  levels)  and  the  possibility  of  discriminating  against  the  scattered 
component  of  the  energy.  The  model  is  stochastic,  producing  a  statistical 
result  baaed  on  the  statistics  of  the  bathymetry.  It  is  also  three- 
dimensional  as  it  incorporates  scattering  out  of  a  vertical  plane.  The 
acoustic  field  is  propagated  from  the  source  to  the  vicinity  of  the 
boundary  interaction  using  a  split-step  parabolic-equation  calculation 
where  the  coherence  function  is  formed.  The  coherence  function  is  deter¬ 
mined  at  the  array  by  propagating  along  characteristic  curves.  A  modified 
Kirchhoff  formulation  is  used  in  weighting  the  scattering  directions  at 
the  points  of  bottom  interaction.  We  have  developed  a  fast  and  accurate 
way  to  estimate  the  large  number  of  scattering  integrals  which  must  be 
evaluated. 


A  flflMWI  M. 


p:  ray  parameter 
X(p) :  offset  function 


LINEAR  GRADIENTS  In  C(z) 

dX(p)  _  2,  x  1  1 

'  P  layers  ^  cos(6i,k 


bk:  gradient  in  layer  k 
bt:  gradient  in  turning  point 

ei,k:  an9le  at  top  layer  k 

ebyk:  an3le  at  bottoni  lQyer  k 
0f  angle  at  top  of  turning  laver 
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VLF  PROPAGATION  #2 

ON  THE  RELATIVE  AMPLITUDES  BETWEEN  PRIMARY  AND  MULTIPLE 
SIGNALS  FROM  SEISMIC  REFRACTIONS  IN  OCEANIC  CRUST 


In  deep  ocean  seismic  re-fraction  experiments  the  energy  in 
the  primary  signal  is  o-ften  less  than  that  in  the  multiple 
arriving  at  approximatel y  twice  the  travel  time.  This  has 
perplexed  interpreters  since  the  multiple  has  a  path  length 
greater  than  the  primary  and  its  amplitude  presumably  should  be 
less.  High  resolution  si  owness/ travel -time  analysis  o-f  an 
expanding  spread  data  set  suggests  an  explanation  in  terms  o-f 
differences  in  the  sound  speed  gradients  in  the  oceanic  crust. 
The  gradient  in  the  upper  part  of  the  crust  is  typically  much 
greater  than  that  of  the  deeper  crust;  therefore,  energy 
refracting  in  the  upper  crust  is  focused  more  and  has  a  decreased 
geometrical  offset.  The  multiple  then  should  have  a  lower  phase 
velocity  and  a  higher  amplitude  over  a  range  window  when  it  is 
refracting  in  the  higher  gradient  part  of  the  crust.  Synthetic 
modeling  using  the  WKBJ  seismogram  method  and  the  amplitude  and 
velocity  analysis  of  the  ESP  data  set  support  this  hypothesis. 


OCEAN  REFRACTION  AND 
BATHYMETRIC  SCATTERING  MODEL  (ORBS) 
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PROPAGATION  BY  COMBINATION  OF  SPLIT-STEP  PARABOUC  EQUATION 
METHOD  AND  COHERENCE  FUNCTION  ALONG  CHARACTERISTICS. 


EXPONENTIAL  SUBSTITUTION 


EXPONENTIAL  SUBSTITUTION  REDUCES  EXECUTION  TIME  BETWEEN  ONE 
AND  TWO  ORDERS  OF  MAGNITUDE  WITHOUT  SIGNIFICANT  LOSS  OF 
ACCURACY. 


EXPONENTIAL  SUBSTITUTION  APPROXIMATION 


EXPONENTIAL  SUBSTITUTION  APPROXIMATION  (COMTO) 


v  '•  V 


A<*  M 


DIFFERENTIAL  CROSS  SECTION 

KIRCHHOff 

X  ”  50.0  m  <7  ”  20  m 

cpto  "  20  deg  L  **  400  m 

W  —  (R/UK^R/L)  HARO  SURFACE 

S  (  20.0,  0.0)  -  28.48 
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DIFFERENTIAL  CROSS  SECTION 


X  “  50.0  m  a  =■  20  m 

< Pi„  ="  20  deg  L  =“  400  m 

W  *  (R/L)K,(R/L)  HARD  SURFACE 


2  (  20.0.  0.0)  -  28.95 


BOTTOM  8LOPE  MOUNTED  ARRAY  ENVIR 
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AD-A154  652  VLF  (VERV  LOU  FREQUENCV)  UORKSHOP:  HELD  AT  SRN  DIEGO 
CALIFORNIA  ON  24-25..  (U)  SCRIPPS  INSTITUTION  OF 
OCEANOGRAPHV  LA  JOLLA  CA  MARINE  PHVSIC.  .  U  S  HODGKISS 
UNCLASSIFIED  1985  HPL-U-11/85  N08014-80-C-0228  F/G  5/1  NL 
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MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- I963-A 


Horizontal  Array  (Rough  Slope) 
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Horizontal  Array  (Rough  Slope) 
Long  Range 


Vertical  Array  (Rough  Slope) 


response  =  0  <JB 


Array  Signal  Gain  of  24-Wavelength  Array  vs. 
Roughness  and  Correlation  Length  of  Isotropic  Bottom 


Normalized  rms  Roughness  a! A 


Beam  Width  of  24-Wavelength  Array  vs. 
Roughness  and  Correlation  Length  of  Isotropic  Bottom 


FUTURE  DIRECTIONS 


■LING  SUPPORT  FOR  CAPE  FEAR  VLF  EXERCISE 


W.  TUNNELL  AND  H.  F.  SCHREINER  .  JR. 
tVAL  OCEAN  RESEARCH  AND  DEVELOPMENT  ACTIVITY  .  CODE  240 
STL  .  MS  39529 

>•  Cape  Fear  Very  Lou  Frequency  <VLF)  Exercise  is  scheduled  off  the 
it  at  Cape  Fear  .  NC  during  June  4  to  July  1  »  1985  as  part  of  the 
5A  6.1  Very  Lou  Frequency  Acoustic  Propagation  Special  Focus  Program, 
exercise  is  a  cooperative  effort  betueen  the  United  States  Geological 
/ey  at  Woods  Hole  and  NORDA.  The  USNS  Lynch^one  other  ship  uill 
port  the  exercise.  Data  acquisition  systems  uill  consist  of  12  ocean 
tom  seismometers  and  a  vertical  hydrophone  array  (  15  element  .  20  m 
cing  )  built  uith  the  VEKA-II  technology.  Underuater  explosives  (  *'*100 
/shot  )  and  a  Mark  6b  sound  source  (  187  dB  at  10  Hz  )  uill  be  used  as 
rces. 

he  modeling  results  presented  here  are  the  array  response  predictions 
culated  uith  the  maximum  entropy  method  using  the  acoustic  field 
dieted  by  the  parabolic  equation  model  solved  uith  the  implicit  finite 
ference  techniques  (  IFD  ).  The  sound  speed  profile  along  a  previous 
5  track  (  USGS  MCL  32)  and  historical  uater  column  data  are  used, 
ilton's  results  <  E. L.  Hamilton,  "Geoacoustic  modeling  of  the  sea  floor" 
.  Acoust.  Soc.  Am.  68(5).  pp.  1313  -1340  »  Nov.  1880  )  for  density  . 
pressional  attenuation  .  shear  sound  speed  .  and  shear  attenuation 
ues  are  used  as  input.  These  inputs  are  used  uith  Kutschale's  version 
the  fast  field  program  to  derive  an  effective  compressional  attenuation 
t  may  be  used  as  input  to  the  IFD  (  the  IFD  does  not  describe  shear 
es  in  the  subbottom).  This  procedure  assumes  that  the  range  dependence 
the  enviroment  has  a  greater  effect  on  the  propagation  than  shear  uaves 
that  shear  effects  can  be  approximated  by  an  increase  in  the  value  of 
pressional  attenuation.  IFD  intensity  contours  are  plotted  for  ranges 
100  to  175  km  (  left  to  right  ).  The  dark  line  indicates  the  subottom 
ch  varies  from  1800  m  deep  to  370  m  deep  left  to  right.  A  high  amount 

mode  drop  out  is  shown  over  the  100  to  140  km  ranges.  The  IFD  results 

used  in  a  maxium  entropy  beamformer  and  the  3  dB  contour  intervals 
obtained  for  arrival  angles  (  -90  to  90  )  at  ranges  from  (A)  100  to 

l  km  .  (B)  80  to  155  km  ,  (C)  60  to  135  km  .  and  (D)  40  to  115  km  . 

se  four  results  cover  the  same  portion  of  USGS  MCL  32  as  the  source 
moved  upslopr  The  results  show  that  the  energy  in  the  mode  drop  out 
ion  is  large  »  incoherent  and  that  only  two  modes  survive  the  modp 
0  out  region. 
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Near  Term  SSQ— 53B  Improveuaeiit 


ABSTRACT 


Summary  of  VLF  Sonobuoy  Program  at  the 


Naval  Air  Development  Center 

James  F.  McEachern 
Code  3043 

Naval  Air  Development  Center 
Warminster,  PA  18974-5000 

NADC's  efforts  in  Very  Low  Frequency  (VLF)  sonar  evolved  out  of  on  going 
improvement  programs  in  passive  sonobuoys.  Dual  isolation  sonobuoy  suspensions 
and  research  in  cable  strunming  attenuation  advanced  the  state  of  the  art  in 
sonobuoy  self  noise  reduction  in  the  areas  of  surface  wave  -induced  noise  and 
flow  induced  vibration  noise.  A  program  of  low  velocity  flow  noise  studies  on 
gradient  hydrophones  led  to  the  development  of  an  effecive,  packageable  flow 
shield  for  gradient  hydrophones. 

The  Infrasonic  Investigation  program  for  target  detection  and  classifica¬ 
tion  drew  on  NADC's  low  frequency  sonobuoy  technology  to  develop  a  sonobuoy  sen¬ 
sor  capable  of  providing  useful  submarine  detection  and  classification  in  the 
VLF  band.  Under  this  program  omnidirectional  and  directional  VLF  sonobuoys 
were  used  to  gather  target  data,  ambient  noise  data  and  to  evaluate  VLF  propa¬ 
gation  models.  The  impact  of  VLF  technology  on  all  airborne  processors  and  on 
sonobuoy  production  and  testing  was  assessed.  A  recommended  self  noise  limit 
for  VLF  sonobuoys  was  generated. 


MODELING  PROCEDURE 


*  USGS  MCL  32  :  COMPRESS.  SOUND  SPEED 


*  HAMILTON  :  DENSITY 

:  COMPRESSION AL  ATTEN. 

:  SHEAR  SOUND  SPEED 
:  SHEAR  ATTENUATION 

* 

*  FFP:  EFFECTIVE  COMPRESSIONAL  ATTEN. 

*  USE  EFFECTIVE  COMPRESSIONAL  ATTEN., 
DENSITY.  AND  COMPRESSIONAL  SOUND 
SPEED  AS  INPUT  TO  PARABOLIC  EQ. 


MAXIUM  ENTROPY  PARAMETERS 


*  INPUT  :  15  ELEMENT  ARRAY 

:  20  M 

:  40  M  ARRAY  DEPTH 

*  MAXIUM  ENTROPY  PARAMETERS 

-  1 0  TH  ORDER 

-  100  AVERAGES  (  25  M  ) 
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FIRST  FIELD  EXERCISE 
JUNE  4  -  JULY  1,  1985  CAPE  FEAR 


PARTICIPANTS:  USGS  AT  WOODS  HOLE 

NORDA  (  SCRIPPS  .  NRL  ) 

ASSETS:  USNS  LYNCH  (  NORDA  ) 

CONTRACT  SHIP  (  NORDA  ) 

12  OBS  UNITS  (  USGS  .  NORDA  ) 
VEKA  (  NORDA  ) 

EXPLOSIVES  (  NORDA  ) 

CW  (  NRL  ,  NORDA  FUNDS  ) 
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TRANSMISSION  LOSS  TO  A  SEISMOMETER 
IN  OCEANIC  BASEMENT 


R,  F.  HENRICK/  J.  R.  ROTTIER 


THE  JOHNS  HOPKINS  UNIVERSITY 
APPLIED  PHYSICS  LABORATORY 
LAUREL/  MD 


TRANSMISSION  LOSS  TO  SEISMOMETER 
IN  OCEANIC  BASEMENT 


#  PNA  (T)  24,470  *  1/18/85 


1C  PROFILE 


FFP  MODEL 


#  PNA  (T)  24,469  *  1/18/85 


"J"  J.I 


J 


J 

D 

C 


N 


LO 

oi 

PM 


n 

3 

3 

J 

X 

3 

u 

D 


n 

u 


m 

E 

u 

"a 


ffi 

3 

Ul 

a 


d  PH  P)  (O  PI  PI  PI  d  P^  PI  P^  P>  PI  PI  p) 

oooooooooooooooooo 


d  m 

O  P^  \6  vfi  i/>  tA 


00  «  O'  O'  ON 


<MNNCM<SNNN<SN 


NN  NN 

mnnnnnoo  oo 

oooooooooo 

C'O'^'OOOOOOOOOO 

ooooooooooooo 


II 

'c 

— 4 

—4 

•4 

—4 

c 

0 

0  0 

0 

O 

0 

0 

O 

0 

0 

0 

O  O 

O 

0 

0  0 

0 

0 

0 

0  0 

O 

O 

0 

0 

O 

0 

O 

O 

O  O 

0 

0 

O  O 

0 

0 

CM  CM  CM  cm 

•  *,•  • 

Go 

0 

0  0 

0 

O 

0 

0 

0 

0 

0 

O 

O  O 

O 

0 

O  O 

0 

0  <*  <t 

-T  CM  CM  CM  CM  CM  CM  O  O  O  O 

*o 

0 

0  0 

O 

O 

0 

0 

0 

0 

0 

0 

O  O 

0 

0 

0  0 

0 

Odd 

«~>oooooooooo 

0 

0  0 

0 

O 

O 

0 

0 

O 

0 

0 

O  O 

O 

0 

O  O 

0 

poo 

00000000000 

u 

a 

9 

0 

• 

0  0 

•  • 

O 

• 

O 

• 

O 

• 

0 

• 

O 

• 

O 

• 

O 

• 

0 

• 

O  O 

•  • 

O 

• 

O 

• 

O  O 

•  • 

0 

• 

000 

•  •  • 

00000000000 

z 

•4^ 

LU 

-S5 

0  0 

•  • 

oomoooooooo 
•  •••••••••• 

V’/-“ 

E 

0  0 

PO-J-OOOOOOOO 

«N  O 

0u">*000'»0'»,000 

mJ 

vn 

-4 

■J-J-jniACOOMoONr^ 

«  *,  1 

a 

u 

MNPIPIPIM^PImM 

»  ..’•  J 

LLl 

in 

0 

• 

O  ^ 

•  • 

CM 

• 

vO 

• 

*-• 

• 

CM 

• 

m 

• 

xr 

• 

as 

• 

CM 

• 

CO  d 
•  • 

ao 

• 

sO 

• 

-44  © 
•  • 

ao 

• 

d  O  ao 
•  •  • 

OOrvoOOOOOoO 

•~~S 

fp 

E 

d 

d  no 

0 

d 

ao 

ao 

Mf 

as 

d  <yi 

d 

d  ao 

<*) 

0  0*0 

OO-TOOO  OOOOO 

1 

«■* 

— •  0 

00 

vO 

MO 

'C 

«n 

m 

yO 

sO 

r-  r* 

00 

O' 

0  — 

d 

in  -4  cm 

■  ■ 

in  in 

-a 

>0“  M 

m  in 

m 

in  M0  M3 

CN  PIPIPI^ 

Ub 

u 

^  «4 

«4 

^4 

** 

•4  ^Hl 

— * 

**4 

—4 

^  ^4  ^ 

V'/-' 

a: 

o 

,E 

0 

d  d 

in 

r-» 

0 

0 

rs. 

as 

m 

0  -• 

CM 

as  «n 

CM 

000 

0000000  _•  _• 

L£a 

X 

0 

m  co 

0 

in 

0 

0 

« 

CM 

•a 

r* 

M3 

** 

0  ao 

Cl 

r—  r** 

-»*»-»oooooooo 

J 

L. 

*4  »4 

d 

•n 

MO 

O 

m 

0 

m 

00  nO 

CM 

d 

mt  in 

r** 

M)  >0  ^ 

<N<vi-*»OOOOOOm»/'l 

•-< 

d 

vO  0 

vn 

^4 

rs  M) 

m 

^  in 

ocaascOm^H-nn 

*1 

CL 

^4 

*-4 

CM 

cm  n 

in  m  in 

in>np»rsNGOflO«4i4B.M 

l 

Ul 

3 

#  PNA  (T)  24,474  *  1/18/85 


SS01  NOISSIWSNVHi 


I 


ISSUES 

I  ACCURACY  OF  PROPAGATION  MODEL 
BOREHOLE  VS.  BOTTOM  WATER 
"  ““  VS.  SOUND  CHANNEL 

SHORT  RANGE  VS.  LONG  RANGE 

I  VALIDITY  OF  GEOACOUSTIC  DESCRIPTION 
.  "  SENSITIVITY  TO  VERT  I  CAL  DETAIL 


HORIZONTAL  ISOTROPY 


> 


■SILTS) 

OMPARISOR  OF  MODELED  PROPAGATION  LOSS  AT  IS  Hz  FOR  ROID 
SOLIS  SEDIMENT  DESCRIPTIONS 


•  AT  RATER-SEDIMENT  INTERFACE  ONLY 

VERTICAL  DISPLACEMENT  STRONGLY  AFFECTED 


•  IN  BOREHOLE,  SOLID  SEDIMENT  MODEL  GIVES  IQ  n8 
OR  MORE  ENHANCEMENT,  AND  HORIZONTAL  DISPLACEMENT 
MORE  ATTENUATED  THAN  VERTICAL  DISPLACEMENT 


SENSORS 


DSDP  SITE  581 
M°N  160°E 


1  CLAY 

I _ SEDIMENT 

I  CHERT 

$I3LOnv 

BASALT 


SUBSTRATE 


PROBLEM  DESCRIPTION: 
BOTTOM  IMPEDANCE  BOUNDARY  CONDITION 
FORMULATION  OF  NORMAL  MODE  PROBLEM 


REFLECTION  COEFFICIENT  R(0,U« 
DETERMINES  BOTTOM  BOUNDARY 
CONDITION  FOR  NORMAL  MODE 


OUTLINE 


BOTTOM  IMPEDANCE  BOUNDARY  CONDITION 

NORMAL  MODE  MODEL . 

BOREHOLE  GEOACOUSTICS 
IMPORTANT  PROPAGATION  MECHANISMS 


ISSUES 


VLF  PROPAGATION 

MODELING  TO  AN  . 

OCEAN  BOREHOLE  RECEIVER 

DR.  ROBERT  A.  KOCH 
APPLIED  RESEARCH  LABORATORIES 
UNIVERSITY  OF  TEXAS  AT  AUSTIN 


ONR,  CODE  425  UA 


NORDA,  CODE  110A 


VIF  Propagation  Modeling  to  an  Ocean  Borehole  Receiver 
Robert  A.  Koch 

Applied  Research  Laboratories 
The  University  of  Texas  at  Austin 
Austin,  Texas  76713-6029 

ABSTRACT 

A  normal  mode  analysis,  in  which  the  seafloor  is  represented  by  an 
impedance  boundary  condition  derived  from  the  complex  reflection 
coefficient  for  a  structured  solid.  Is  used  to  predict  the  dominant  features 
of  low  frequency  propagation  from  a  water  column  source  to  near  bottom 
and  subbottom  receivers.  The  propagation  loss  calculations  suggest  that 
propagation  to  a  receiver  In  the  basalt  Is  predominately  via  shear  waves. 
Based  upon  these  and  other  results,  suggestions  for  future  efforts  are 
offered. 
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FREQUENCY 


SITE  581  TL  COMPARISON 
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OBJECTIVE 
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/  APPROACH 

CONDUCT  A  TUOrSEASON  ACOUSTFC/GEOACOUSTIC  EXPERIMENTAL  INVESTIGA 
TION  OF  S/N  IMPROVEMENT  FOR  BORE  ROLE  SENSORS  TN  SHALLOW  MATER* 
SUPPORTED  BY  ACOUSTIC  AND  6EOACOUSTIC  HO DEL INC  AND  GEOPHYSICAL 
MEASUREMENTS. 
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FIGURE  6 

FREQUENCY  DEPENDENCE  OF 
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